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("^ ' Abstract. We present the results from new relativistic hydrodynamic simulations of binary neutron star mergers using realistic 

c7^ , non-zero temperature equations of state. We vary several unknown parameters in the system such as the neutron star (NS) 

CN| ' masses, their spins and the nuclear equation of state. The results are then investigated with special focus on the post-merger 

J, ' torus-remnant system. Observational implications on the Gamma-ray burst (GRB) energetics are discussed and compared 

Oj' with recent observations. 
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INTRODUCTION 

T 1 ■ 

> 

^SJ ' Merger events of binary neutron stars (BNS) and neutron star - black hole (NSh-BH) binaries do not only belong 

^sO ■ to the strongest known sources of gravitational wave (GW) radiation, they are also widely favored as origin of the 

'/^ ' subclass of short, hard Gamma-ray bursts (GRBs) fr, T, ^, M^, ^. The recent first good localizations of short bursts 

^r I by the Swift and Hete satellites were interpreted as a possible confirmation of this hypothesis |6, 7, 8, 9, 10], because 

the bursts have observational characteristics which are different from those of long GRBs, but which are in agreement 
^^ I with expectations from compact object mergers as e.g. the energetics and the association with old stellar populations. 

The central engines of such bursts are still poorly understood and observationally undetermined. But it seems unlikely 
that the energies required for typical short GRBs are set free during the dynamical phase of the merging of two NSs. 
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Instead, the following secular accretion phase in a postmerger system consisting of a central BH and a surrounding 



r system consisting c 
O ' torus seems to be a much more promising source (see e.g., fill flZ \13L UA. UjR ) provided the torus is sufficiently 

massive. Thermal energy release preferentially above the poles of the BH by the annihilation of neutrino-antineutrino 
^ ' (vv) pairs canlead to collimated, highly relativistic jets of baryonic matter 1 161 . The gamma radiation is then produced 

in internal shocks when different blobs of ultrarelativistic matter in the jet collide with each other. When the jet hits 
the ambient interstellar medium, the less energetic afterglow radiation is produced. 

Torus mass as well as BH mass and rotation are thus crucial parameters that determine the energy release from the 
«2j ■ mer ger event. Here, we summarize results from our recent 3D relativistic hydrodamics simulations of BNS mergers 

C^ ' ITtI [TSi] to investigate the dependence of the torus and BH properties on the initial BNS parameters like the nuclear 

EoS, the NS masses and mass ratio and the NS spins. 

In the following, we define the torus by the condition that the specific angular momentum of a fluid element in 
the torus has to be larger than the innermost stable circular orbit of a Kerr black hole with the mass and the rotation 
parameter of the remnant. The latter in turn is defined as all the matter which is not in the torus. 

RESULTS FROM RELATIVISTIC HYDRODYNAMICS SIMULATIONS OF BNS 

MERGERS 

Using our relativisic smooth particle hydrodynamics (SPH) code llSll . we have simulated a large set of BNS merger 
models, varying the EoS, the initial neutron star mass ratio, the total mass and the NS spins. 

We have used two realistic, non-zero temperature EoSs whose supranuclear parts base upon two fundamentally 
different approaches to determine the nuclear interaction. The EoS of Shen et al. II9II20I1 ("Shen-EoS") bases on the 
phenomenological relativistic mean field (RMF) approach, while the EoS of Lattimer and Swesty |21] ("LS-EoS") 
employs the finite-temperature compressible liquid droplet model for the nucleus |22]. In addition, we consider an 
ideal gas EoS whose adiabatic index F is adjusted to fit the Shen-EoS in the supranuclear regime. 



TABLE 1. Setup parameters and characteristic data of our considered models. Mi and M2 denote thie individual gravita- 
tional masses in isolation while M^um ~ Mi +M2 stands for the sum of the two (Note that the total gravitational mass M is 
slightly smaller than Mg^^m because M also involves the negative gravitational binding energy between the two stars.). My is 
the total baryonic mass, whereas q = Mi /M2 and qM = Mq.i /M0.2 is the gravitational and baryonic mass ratio, respectively. 
'Shen' stands for the full, finite temperature Shen-EoS, 'Ls' denotes the Lattimer-Swesty-EoS and 'ideal gas' means ideal 
gas EoS. The NS spin states are denoted as 'irrot' for irrotating, 'corot' for corotating, 'counter' for counterrotating, 'oppo' 
for oppositely oriented spins and 'tilted' for a spin orientation tilted relative to the orbital spin. More specifically, we choose 
in model S 1414tl fij = and ^2 = 0.041 * (0, 1 , 0) (in cartesian coordinates) which corresponds to a spin period of the sec- 
ond NS of ~ 1 ms, while in model S1414t2, Sli = 0.01 * (0, 1, 1) and SI2 = 0.01 * (0,-1, 1) is chosen. Tmax is the maximum 
temperature in MeV reached in the system during the whole evolution. It is obtained by averaging the particle temperatures 



on a grid of 1.5km side length. 
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JK 



is the (total) spin parameter measured immediately after merger 



and Mtoms denotes the estimated (baryonic) torus mass. Note that in model LS1216 the maximal temperature and the torus 
mass cannot be determined because the merger remnant collapses immediatedy after merging. 
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Initial data is obtained using a relaxation technique which drives the system into an equilibrium velocity field which 
is given by the orbital separation, the NS spins Q.\ and ^2, respectively and the orbital angular velocity il. The latter is 
iteratively determined such that an orbital equilibrium state between gravitational and centrifugal force is established 
(see [18] for an extensive description). The NS spins are chosen either to be zero (irrotating NSs), or to be equal 
(corotating NSs) or opposed (counterrotating NSs) to the orbital spin. A few cases with non-aligned NS spins have 
also been considered (see Tab.^. 

In addition, we assume the initial NSs to be cold and to be in neutrino-less beta-equilibrium which determines the 
initial distribution of the internal energy u and the electron fraction Ye (In the ideal gas case, the polytropic relation 
M = Kp^^^ , with K fitted to the Shen-EoS, is employed to initialize the internal energy). In Tab.ffl we summarize the 
parameters of our models. The simulations are carried out with typically 400'000 SPH particles from slightly outside 
the innermost stable circular orbit, through the late inspiral, merging and torus formation until either the collapse of 
the merger remnant sets in, or a quasi-stationary torus-remnant state has formed. 



In models with a mass ratio considerably different from unity, e.g. in models S107193, S1216 or S1317, the less 
massive but larger star is tidally elongated and disrupted to form a long spiral arm which is largely accreted onto 
its more massive merger partner. This leads to a non-axisymmetric, differentially rotating central object and a thick 
massive torus around it. The core of the central remnant forms out of the more massive NS while the accreted material 
from the less massive partner is wound up to end in the outer shells of the remnant. The material sitting in the spiral 
arm tip receives enough angular momentum due to a 'sling-shot-effect' to escape accretion and to contribute to the 



main part of the future torus. This carries a significant amount of the angular momentum from the remnant out into the 
torus (see Fig.|2j. On the other hand, in case of a mass ratio near unity, as eg. in S1414 or S138142 the two partners 
smoothly plunge together without any disruption and formation of a primary spiral arm. The two original NS cores are 
kept mostly intact during merging and form together a rotating twin-core, bar-like structure at the center of the merger 
remnant. This rotating bar drives the formation of secondary spiral arms at the remnant surface after the dynamical 
merger phase and leads to ejection of material into a small torus. Note that primary spiral arms and large torus masses 
may still form in these systems, if the initial NS spin setup is favourable, e.g. in case of corotating NSs. 
For fixed EoS and NS spin, we find, that the torus mass depends roughly linear on the mass ratio for ^-values larger 
that about 0.8, whereas below that point, the torus mass saturates at about Mtoms — 0.25Mr.). In total, the torus masses 
range from about O.OSM© at g' = 1 to O.25M0 atq — 0.55 (see Fig.^ panel (b)). 

The NS spin setup mainly determines the total amount of angular momentum in the system and thus directly 
influences the large-scale postmerger evolution, the amount of material carried into the torus and the spin parameter 
of the merger remnant. The latter particularly influences the torus mass via our definiton of the torus. If the NS are 
initially corotating, the increase in total angular momentum is about 8% compared to the irrotating case. Similarly, 
counterrotating NSs decrease the total angular momentum by the same amount. 

We find, for fixed EoS (Shen-EoS) and mass ratio {q — 1), that the torus mass is indeed mainly determined by 
the amount of total angular momentum available in the system whereas details in the merger dynamics are of minor 
importance. Models with a different spin setup, but a similar total angular momentum yield comparable torus masses. 
This is the case for models S1414co (corotating) and S1414t2 (Tilted spins summing up to a total spin in z-direction), 
and, to a smaller extent, for models S1414 (no spins) and S1414o (oppositely oriented spins). This suggests that 
the torus mass dependence on the NS spin setup can be approximately reduced to a single parameter, the total angular 
momentum in the system (see Fig.^ panel (a)). In total, we find torus mass values between O.O4M0 for counterrotating 
NSs and 0.3Mq for corotating NSs. 

The EoS influences the merger dynamics and postmerger outcome in different ways. Initially more compact NSs, 
as obtained with a soft or compressible EoS, lead to a longer inspiral phase because the tidal instability sets in at a 
smaller orbital separation. Therefore, more angular momentum is carried away by gravitational wave radiation and 
the angular momentum in the post-merger system is reduced which favours a smaller torus. On the other hand, the 
angular momentum transport through non-axisymmetric motion, as e.g. a rotating bar-like core, or through numerical 
viscosity is found to be more efficient if the remnant is more compact and more rapidly rotating. This in turn favours a 
larger torus. In model LS 1414 (see Fig.^ panel (c)), about twice as much of the angular momentum is radiated away 
(15% compared to 8% in model S 1414) and the torus mass shortly ( 3ms) after merger is about half as large ( O.OSM© 
versus O.O6M0). During the following post-merger evolution the torus grows continously. The ideal gas model P1315 
leads almost to the same torus mass as the Shen-EoS model S 13 15. This suggests that the torus mass is determined by 
the supranuclear part of the EoS. 

IMPLICATIONS FROM OBSERVATIONS OF SHORT GRBS 

In this section we discuss the implications to the postmerger system properties from the recent observations and 
locations of six short GRBs by the Swift and HETE satellites. Several observational results suggest, that at least part 
of the short GRBs may come from binary NS mergers. Short GRBs seem to be partially associated with elliptical 
galaxies where star formation has stopped long ago. They also seem not to be associated with supernovae. This is 
consistent with the long inspiral evolution of binary NS systems prior to merger (see e.g. |23]). 

From the measured isotropic equivalent energy Eyjso, we can now estimate, within the merger model, the mass 
which is accreted onto the central BH: 

Ey.iso — /l/2/3/4/n lA/accreted 

Here, /i is the efficiency at which the accreted rest mass energy is converted to neutrino emission, /2 is the conversion 
efficiency of neutrino-antineutrino annilihation to e* pairs, /s is the fraction of the e*-photon fireball energy which 
drives the ultrarelativistic outflow with Lorentz factors F > 100 as required by GRBs, /4 is the fraction of the energy 
in ultrarelativistic jet matter which can be converted to gamma radiation in dissipative processes that occur in internal 
shocks in the jet and fa = 2iijet/(4;r) = 1 — cos 0jet denotes the jet collimation factor defined as the fraction of the sky 
covered by the two polar jets (with semi-opening angles 0jet and solid angles iijet). 
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FIGURE 1. (a) Torus masses versus total spin parameter flmerger = ■'merger /A/mereer °f ^^^ system shortly after merger. Plotted 
are two series of models, with fixed EoS and with mass ratio ^=1 ('+'-signs) and q=0J5 ('x'-signs), respectively.(b) Torus masses 
versus mass ratio q. Three different values of the total gravitational mass are chosen and marked with 'x'-sigs (M^^m = 2.6M0), 
circles (Msum = 2.8M,3) and '+'-signs (Msum = 3.0M(.,). These models use the Shen-EoS and are initially irrotating. (c) Torus 
masses for different EoSs. The Shen-EoS models are plotted in red, while the LS-EoS models are colored in blue and the ideal gas 
EoS model in black. The more compact LS-EoS model LS1414 leads to a significantly smaller torus mass then the corresponding 
Shen-EoS model S1414 shortly after merger. During the subsequent evolution, angular momentum transport by non-axisymmetric 
tidal interaction and numerical viscosity increases the torus mass. Note that the peaks initially after merger are an artifact of the 
non-linear feedback between torus mass, remnant mass and the angular momentum of the last stable orbit in our torus-determination 
algorithm. 



TABLE 2. Estimated accreted masses, M^^c, and lower bounds for 
the average mass accretion rates. The quantity £y,iso is the isotropic- 
equivalent y-ray burst energy (corrected for the cosmological redshift 
z of the burst), and ty is the GRB duration at the source, computed 
as ty = 790/(1 + z) from the measured 90%-inclusive interval of high 
energy emission. The observational data were taken from Fox et al. 
and Soderberg et al. [36.1 . 
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"Typical" values from merger and accretion simulations are: /i - 0.05 1 15, 24], /2 -■ 0.001 ... 0.01 1 12, 24*], /s -- 0. 1, 
/n ~ 0.01 ...0.05 L16.1 and /4 < 0.2 L25>.26>.27> and references therein]. To estimate the accreted mass, we use the set 



(/i,/2,/3,/a,/4)- (0.1,0.01,0.1,0.01,0.1) 



(1) 



to obtain the values in Tab. |2l In addition, we can calculate the average accretion rate Mace. = ^acc/facc ^ M^cc/ty 
assuming that the observed GRB duration ty is an upper limit to the accretion time face fl^- Despite of the crude 
assumption in Q, we find in nearly all cases values for the torus mass in the range obtained from our theoretical 
models. In addition, the estimated mass accretion rates lie in the ballpark of the results from postmerger accretion 
disk simulations LIZ .28- .24- ■15.1 . Note however that all the f- factors have large uncertainties and are likely not to be 
constant but to depend on the torus mass. For example, the neutrino annilihation efficiencies sensitively depend on the 
temperature which in turn depends on the fluid density and thus on the torus mass. Therefore, the large values for the 
torus mass and the accretion rate from the very luminous burst GRB051221 are not astonishing but suggest different 
values for the f-factors or, alternatively, a NSh-BH merger as a source. Newtonian simulations of NSh-BH mergers 
ll29ll30tl3lll32Ll33ll34ll35ll suggest that, for certain mass ratios andEoSs, large torus masses in the range of 0. 1 ... IM© 
may be formed. 
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FIGURE 2. Characteristic snapsiiots of four representative models. Plotted is every 10th SPH particle in a slice around the 
equatorial plane. The matter of the two NSs is color-coded by green and blue particles. Particles ending up in the torus are plotted 
in red while the yellow circles represent particles which currently have an angular momentum larger than the angular momentum 
of the last stable orbit of the merger remnant at the end of the simulation. (By definition, the yellow and red particles coincide at 
the end of the run.). Panel (a) shows the development of secondary spiral arms in model S1414. Panel (b) shows the large, primary 
spiral arms in model S1414co. Panel (c) shows the post-merger evolution of model L1414. The merger remnant in much more 
compact and no secondary spiral arms are forming. As a consequence, the torus is much smaller and is developing on a longer 
timescale. In panel (d), we see the formation of the large primary spiral arm in model S1216 
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